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In this work, we solve exactly a Kondo lattice model in the thermodynamic limit. 
The system consists of an electronic conduction band described by unconstrained 
hopping matrix elements between the lattice sites. The conducting electrons interact 
with a localized impurity spin at each lattice cell. We have found the exact thermody- 
namics, the ground state energies of the system. At T = 0, we explicitly demonstrate 
that the system exhibits a metal-insulator phase transition at half-filling. In the limit 
of strong coupling between the impurity spin and the conduction electrons, J = +00, 
we have solved the system on a lattice of any size L. The ground states are the RVB- 
type Jastrow product wavefunctions. Various correlation functions can be computed 
for the impurity spins, and for the singlets formed by electrons and impurities. 
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The heavy-fermion systems have been of great theoretical interests in recent years |T] || . 
A basic model for a heavy fermion system is an S = 1/2 Kondo lattice with a conducting 
band interacting with a localized impurity spin at each unit cell. Quite recently, a possi- 
bility has been discussed that heavy-fermion superconductors involve odd-frequency triplet 
pairing, through a Majorana representation for the local moments to avoid a Gutzwiller pro- 
jection ||. In spite of several numerical study and approximate approaches to the Kondo 
lattice model, it is too difficult to find the exact solutions for the system. 

In one dimension, the Kondo lattice model has been studied through various numerical 
and approximate methods. The results indicate the existence of a finite spin gap and a 
finite charge gap simultaneously at half filling [f^-|ll|. The numerical renormalization group 



method can even deal with a finite lattice up to 24 sites |10[ , confirming further the charge 
gap at half-filling in the thermodynamic limit. It also shows that the charge gap is larger 
than the spin gap at half-filling P|§. The finite charge gap shows that the system described 
by the Kondo lattice model at half-filling is an insulator. However, due to the complexity, it 
has been impossible to prove the existence of the metal-insulator phase transition at the half- 
filling through any exact approach without employing an approximate scheme or numerical 
methods. 

The single impurity Kondo model can be reduced to a one dimensional problem, as only 
the S-wave scattering off the impurity spin is important. The single impurity Hamiltonian, 
which represents a local Fermi liquid, is exactly solvable with the Bethe ansatz if the con- 
duction band is described by continuum electrons with linear spectrum at the Fermi surface 



T3|]l4}| . The Bethe-ansatz solutions explicitly demonstrate the dynamical energy scale that 



enters the expression of the thermodynamics of the system [|12|-|T4]1 . The screening of the im- 
purity spin by conduction electrons, the crossover from weak coupling fixed point to strong 
coupling fixed point behavior can be explicitly illustrated with the Bethe- ansatz solutions 

BE 

Despite the Bethe-ansatz solutions for the single impurity Kondo model, no exact solution 
for more than one impurity spin is available by far, either in continuum case or on a lattice. In 



this work, we introduce a Kondo lattice model in which the conduction electrons hop with 
unconstrained hopping matrix elements. The local moment at each site interacting with 
the conducting electrons is described by a spin 1/2. We solve this Kondo lattice system 
exactly in the thermodynamic limit at any temperature. Although the hopping matrix of 
the conducting electrons is far from the realistic case ( electrons hop between the nearest 
neighboring sites in the tight binding picture), our model provides the first example where 
its exact solutions explicitly demonstrate the metal-insulator phase transition at half-filling, 
as the interaction between the local impurity moments and the conduction band is turned 
on. 

The Hamiltonian of the system is given by: 

— , — , L 8 

H = Z *w4rC*«r + JY, cta-lfciP ' (*)) W 

l<i^j<L(T=1,i i=l 

where the hopping matrix element for the conduction electrons is given = —t for all 
i j£ j. J is the coupling constant between the local impurity moments and the conduct- 
ing electrons. We assume that t > 0. The electron creation and annihilation operators 
have the usual anticommutation relations {c ia , Cjp} + = 0,{c\ a ,Cjp} + = &ij5 a f3- The local 
moments are described by the spin 1/2 operators, that is, [Sf(k),SH(k)] = iSf(k) (plus 
two other commutation relations obtained by the cyclic permutations of x,y,z), with the 
relation S"j(k) =3/4, for all the sites k = 1,2, • ■ ■ ,L. Because of the special form of the 
conduction band, the dimensionality of the lattice is irrelevant, and the system is basically 
one dimensional. In the following, we always discuss the thermodynamic limit where the 
lattice size L — ► oo. 

The Hamiltonian of the system has 577(2) spin symmetry and the 577(2) isospin sym- 
metry, if we represent the local impurity moment through an / fermion: Sf(k) = ft a -^-fk/3 
at each site k with the condition that J2 a flafka = 1- The isospin operators are given by 

i + = Ei(-iY(4Ai - 44)' r = ( Jt ) f > J * = £*(cU* - 4/- - !)/ 2 0- In the Fourier 

space, the conduction electron Hamiltonian takes the following form 
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tij C icr C . 

l<i^j<La=t,l 



J2 £ e (*O c L c *. 



(2) 



where the Fourier transform of the fermion operators are c ( 



,t 



_J_\^L ikx ir \ 

JJ72 Z^j=l c <-"i<T> L fecr 



7772 Y^i=i e~ ikXi Cj ta . The wave vectors are given by k = 2nir/L, with n = — (L — 1)/2, ■ ■ ■ , (L — 
l)/2 . The dispersion of the conduction band is e(k) = (—t'L)5kfl + 1- Here, £' = £, how- 
ever, we treat this as an independent parameter in the following. The dispersion indicates 
that there are energy levels (—t'L + 1) and t, and the plane waves with energy t have large 
degeneracies. Because of this, to find the thermodynamics of this Kondo lattice model, we 
shall employ the argument of van Dongen and Vollhardt [[RJ for the Hubbard model. The 
impurity spin interaction Hf m = JY^=\Aor^t c ifi ' Sf(i) can be treated as a perturbation, 
and the thermodynamic grand potential may be expanded as a power series of the coupling 
constant J. In the thermodynamic limit, we may see that the grand thermodynamic poten- 
tial comes from two independent parts, with a correction of relative order 0(1/ L) or 0(1), 
which is ignorable in the thermodynamic limit. 
The partition function is given by 



where N the number of electrons, fi the chemical potential, and the grand thermodynamic 
potential is Q = — In Z. The trace is over the Hilbert space for electrons and im- 
purities. Treating Hf m as a perturbation term, we write Q = Q — /3 _1 Wt'(/9, J), where 
Qq = — f3~ x lnTrexp(— f3(H c — fiN)). The contributions from all the connected diagrams 
are given by 



Here, the notation of the expectation value < TA(ti)A(t2) ■ ■ ■ A(r n ) >q is defined to be the 
sum of all the connected diagrams of 



< TA(t,) ■ ■■A(T n ) >= TrTe~ S« dT(Hc -» N) A^A^i) ■ ■ ■ A(r n ) /Tre~ ^ d < H ^ N \ (5) 



The leading contributions to the quantity W#(/3,J) only come from those connected dia- 
grams consisting of the bare electron propagators with momenta k ^ scattering off the 
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(3) 




(4) 



impurity spins. The electron propagator with momentum k = has t'-dependence as Lt', 
so that its contribution to the perturbation expansion is of order of 0(1/L) or 0(1). The 
propagator with momentum k 7^ is independent of the momentum, as well as independent 
of the parameter t'. Therefore, we can put t' — in the Eq. (^), and then include the 
momentum k = in any of the momentum integral. The second step only gives rise to an 
error of order 1, so that 



Wt> = Wv = o + 0{l), (6) 

where the first term is of the order 0(L). Therefore we see that in the thermodynamic limit, 
Wt> is independent of the parameter t'. 

According to the above argument, the grand thermodynamic potential comes from two 
contributions that are completely decoupled. Both parts can be found in a straightforward 
way, and we have in the thermodynamic limit 

u = lim ^ = -2t - /T 1 ln[2 + z(e 3/3J / 4 + 3e~^ 4 ) + 2z 2 ] (7) 

L — >oo J_j 

where z = and \i = \i — t. The electron density n of the infinite system, as function of 
the chemical potential \x is given by n = — (d(j/dfj,)p, i.e., 

U ~ ' [2 + z(e 3 M* + 3e-? J / 4 ) + 2z*] ' (8) 
Given electron density n, the renormalized chemical potential p, is found from 

e M = 1 r_ (1 _ n)(e W4 + 3e -/JJ/4) + 

+ [(1 - n) 2 (e 3/3J/4 + Ze- f3J l A ) 2 + 8n(4 - 2n)f 2 }. (9) 



At half- filling n = 1, and thus /i = 0. 

a/3 im- 



Let us then calculate the energy density for the infinite system, i.e., e = jin + 



we obtain 



o ^ ^[(3J/4 + /l)e^/ 4 + 3(-J/4 + / 2)e-^/ 4 ]+4^ 2 1 



For fixed number of electron density n, taking the low temperature limit (3 — » oo, the 
above becomes the ground state energy of the system. Let us first consider the case where 
the interaction between the local moment and conduction band is antiferromagnetic, J > 0. 
At less than half-filling n < 1, we see that fx — > —3 J/4 when /3 — ► oo. The ground state 
energy is thus 

3 

e G = --Jn-2t + tn. (11) 

At half-filling n — 1, p, — 0, and the ground state energy is 

e G = --J-2t + t. (12) 

At more than half-filling n > 1, the chemical potential /2 — > 3 J/4 when /5 — > oo, and the 
ground state energy is given by 

3 3 
e G = --Jn-2t + tn + -J(n-1). (13) 

From these results, we thus obtain 

dec 3 J 

dec . 3 J 

^ Ul+ = T + t 

dec 3 J 

^Un = - T +t,n <l 

dec i 3 J . s 

-l—\n-^n = -r + t,n >l. (14) 
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To study the nature of the ground state, we employ the idea due to Mattis |L7j|18]], defining 
ll + = E G (N+l)-E G (N),fi- = E G (N)-E G (N-1). The idea is that the system is a metal 
if 5li = li + — \l~ = 0, while the system is an insulator if 8/i > 0. At half-filling, there is a 
kink in the chemical potential ll + > ll~ , indicating that the system is an insulator due to the 
impurity spins. At less than half-filling n < 1, jl + = ilr = —3 J/4, and the system is metallic. 
When the filling number is larger than one, n > 1, ji + = ilr = 3 J/4, indicating that the 
system is metallic. From these analysis, we see that a metal-insulator phase transition occurs 
at half-filling, as we turn on the antiferromagnetic interaction between the local impurity 
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moments and the conduction electrons. This phase transition would occur even at very- 
small interaction parameter. 

Finally, the analysis of the situation of ferromagnetic coupling J < leads to the same 
conclusion: At less than half-filling n < 1, p, — > — | J|/4 when (3 — > oo, and the ground state 
energy density is 

e G = -\J\ n /A + {-2t} + tn. (15) 

At exact half-filling, the chemical potential Jx — 0, and the ground state energy of the system 
is given by 

e G = {-2t-\J\/A} + t. (16) 

When the electron density is more than 1, we find that ft, — > |J|/4 when /3 — > oo, and the 
ground state energy is 

e G = \J\n/4 + {-2t- \J\/2} + tn. (17) 

From these results, we see that the system is metallic if n ^ 1, as indicated by 8ji = 
fi + — fx~ = 0. However, at half-filling n — 1, we have 5/x = |J|/2, showing that the system 
is an insulator. 

At zero temperature, when there is no interaction J = 0, the system is metallic ( a 
simple Fermi liquid ). This fixed point is unstable at half-filling against infinitesimal small 
interaction. Our analysis has shown that, at half-filling, a finite charge gap of the order 
| J | is open in the excitation spectrum, when turning on the interaction. As we have seen, 
due to the long range aspect of the hopping matrix, for the thermal potential, the kinetic 
energy part decouples with the interaction energy in the thermodynamic limit. Loosely 
speaking, at half-filling, to lower the energy of the system, each conduction electron would 
attempt to form a singlet ( triplet ) at each site, when an antiferromagnetic ( ferromagnetic 
) interaction is turned on. For an electron to hop from a site to any other site on the chain, 
even with the long range hopping, it would have to break two singlet ( triplet ) pairs first, 
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giving rise to a finite charge gap proportional to \J\. In the end, both ferromagnetic and 
antiferromagnetic interactions will drive the system to an insulating phase. In the above 
study, the diagrammatic argument does not give us any information on whether one can solve 
the Kondo lattice model on a finite size lattice, nor does it provide explicit wavefunctions 
of the system. However, as will be demonstrated below, in the limit of strong interaction, 
one can actually find the ground state wavefunctions rigorously, for a lattice of any size L. 

In the case t < 0, at zero temperature, similar analysis indicates the metal-insulator phase 
transition occurs in the system at half-filling, when one turns on J. At finite temperature 
T, the free energy density is found to be 



where z = e^~^ and /i is the chemical potential. The diagrammatic analysis has shown 
that the free energy of the system can be written in closed form. However, we would like to 
emphases that one can not treat the system as if the hopping matrix were t = 0, and each 
electron would form localized singlet ( triplet) with the impurity at each site for positive 
J ( or negative J), uncorrelated with each other. The long range hopping matrix would 
delocalize them, and induce strong correlations between them, as shown below. Consider 
the special situation t < 0, J = +oo. In the limit J = +oo, supposing that there are N e 
electrons on the lattice L, with N e < L, each electron will attempt to form a singlet with 
the impurity spin at each site, to lower the energy of the system as much as possible. Some 
unpaired impurity spins are left over on the lattice. Therefore, we work in the Hilbert space 
where each site can be either a unpaired impurity spin or a singlet of electron-impurity 
bound state. Due to the hopping of the conduction electrons, the singlets can hop on the 
lattice. Let us rewrite the Hamiltonian in the following way: 




(18) 




(19) 



In the limit J 



+oo, the basis vectors can be written as 



a > 



2 ~N e /2 TTQ-P a ) C f C f •••C 

Z 111 1 r nfc) C X l7l C X 2 72 U - 



> 



Li=l 



8 



<S> °"2, • • • , Pi, • • ■ , fo, ■ • • , &L-N e >, (20) 

where the singlets are located at positions {x} — (x± < x-i < • ■ • < xjyj, the unpaired 
impurity spins (a±, 02, • ■ • , <JL-N e ) are positioned at sites {y} — (y± < y-i < ■ • • < yL-N e )- 
Here, the operator P 7i ^ permutes the spin indices 7« and to form a singlet of electron 
and impurity at site Xi. 

Denoting the projection operation onto the J = +00 subspace by P, the projected 
Hamiltonian would take the form: 

H = PHP = PTP + c = H l + c ) (21) 

where T is the kinetic energy of the conduction electrons, the infinite constant c = (—J/4)N e 
only shifts the origin of the energy of the system, a reference energy which is unimportant 
physically. In the space where the z-component of the total spin is fixed, that is, S z = M, 
the number of the unpaired up-spin impurities is A = M + (L — N e )/2, the number of the 
unpaired down-spin impurities is B = —M+(L — N e )/2. The dimension of the Hilbert space 
associated with J = +00 is thus C^ e x Cf i _ Ne . Any eigenenergy state of the Hamiltonian 
Hi = PTP can be written as a linear combination of the basis vectors, 

\<j>>=J2C(a)\a> . (22) 

a 

To find the ground state wavef unctions, we may identify the singlets as spinless fermions, 
the unpaired impurities as hard core spin 1/2 bosons hopping on the lattice. Let us consider 
a system described by the following Hamiltonian: 

h = (1/2) £ P G (t ijg l gj b ia bl)P G (23) 

where the b fields are bosonic, g fields are fermionic, and tij = \t\ for any i 7^ j. The b fields 
commute with the g fields, and J2i9i9i — N e . The Gutzwiller projector P G restricts the 
system to be in the subspace g\gi + J2a=i,l &]<A<t = 1 for i = 1, 2, • • • , L. For this system, 
the basis vectors may be represented as follows: 
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I« >= <4<4 ■ ■ ■ 9l N A<n b i<r* ■ ■ ■ b L- N ^ L . Ne 10 > • (24) 
One can easily verify the following matrix elements 

</?|#i|a>=</3|/i|a>. (25) 

Therefore, the two systems considered here are isomorphic to each other, and we have 
the one-to-one correspondence \a >^ \a > for the basis vectors. Using the superalgebra 
representation 

P F (i)FlP F (i) = P G bl 9l P G 

P F (i)F la P F (i) = P G b ia gjP G , (26) 

where P F (i) = (1 — FtF^F^Fii) is the Gutzwiller operator, and F is a fermionic field, we 
have 

h = (-|t|/2) J2 P F FlF ja P F , (27) 
No- 
where P F = Yif-i P F {i), an d the number of the F fermions on the chain is N F = L — N e . 
For the system h = (— 1£|/2) Ystfja PFF} a Fj a P F , it has been proved before that the 



ground state energy of the Hamiltonian with N F = L — N e is [[0|] 



e G = (-|*|)iV e . (28) 
This ground state is highly degenerate, and one has 

L-N e -2 

\<fo >= PfF^I II ^Ll°>> ( 29 ) 

i=l,fc i ^0 

where F^pF^ and Fl_ a _ creat up spin F fermion with momentum 0, down-spin F fermion 
with momentum 0, and F fermion with spin <7j and momentum respectively One can also 
find some trivial excited energy levels. For example, one eigenenergy is e% = (\t\/2)(L — N e ), 
another one is €2 = (— |t|/2)iV e . 

We can write out some of the ground state wavefunctions explicitly. Denote any state 
vector in the following fashion 
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|0 >= £ <S>({X},{Y})f[F Y ^f[F} [jl F x ^\P>, (30) 
{x},{y} i=i i=i 

where the reference state |P >= n^i-^tl >, <2 = ^e, 5 = — M + (L - iV e )/2, the 
amplitude $ is antisymmetric in the positions {¥}, while symmetric in the positions {X}. 
We have proven that the following Jastrow wavefunctions are the eigenenergy states with 
energy (-|*|)JV e , 

$({X}, {Y}) = ^/mm^Xi+m^Yi) J[^ Xi - Xj) H diYi - Yj) HdiX, - Yj), (31) 

where the function d(n) = sin(7rn/L), and the quantum numbers m s ,mh are integers or 
half-integers, which make sure of the periodic boundary conditions, and satisfy the following 
constrains 

(Q + B + l)/2 <m h <L-(Q + B + l)/2 

(Q + A+ l)/2 < (m h - m s + L/2) < L - (Q + A + l)/2. (32) 

The effect of the up-spin F fermion hopping operator can be calculated readily, when it acts 
on the Jastrow wavefunctions. Rewriting the Jastrow wavefunctions in terms of the positions 
of the up-spin F fermions and the holes, the resultant amplitude still takes a similar Jastrow 
product form. The down-spin F fermion hopping operator effect can then be handled in the 
same way. With these simple calculations, we have been able to show that the wavefunctions 
Eq. (|3~lD are indeed the eigenfunctions of the Hamiltonian h, with the ground state energy 
ec = (—\t\)N e , as long as the quantum numbers m s and m h satisfy the constrains Eq. (|32"|). 

It should be remarked that the singlets and the unpaired impurity spins are strongly 
correlated with each other, as seen from the Jastrow product wavefunctions, due to the long 
range hopping aspect. Various correlations between the unpaired impurities and the singlets 
can be computed exactly in compact form. Given the Gutzwiller- Jastrow wavefunctions 
with the quantum number m s and mh, we can trivially generalize Forrester's work (m s = 
L/2,rrif l = L/2) |2tJ to these wavefunctions. In the conventional case, where the conduction 
electrons hop between the nearest neighboring sites, in the strong interaction limit J = +oo, 
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the general eigenenergies can be found with the Bethe-ansatz solution as for the ID Hubbard 
model with strong repulsion between the F fermions. It was impossible to obtain the closed 
form of the correlation functions of the impurity spins and the singlets exactly, except 
their approximate long distance behaviors using finite size scaling analysis of conformal field 
theory. Finally, we would like to note that these solutions are valid for any lattice size L, 
unlike our first part of finite J discussion that deals with the large lattice. Moreover, like 
in the long range Hubbard model, the thermodynamic limit and the large J limit do not 
commute with each other. 

In summary, we have solved exactly the Kondo lattice model in the thermodynamic limit. 
Due to the long range hopping matrix, we have been able to show that the contribution of 
some scattering processes to the thermal potential scales as 0(1) or 0(1/ L), while the rest, 
which can be computed explicitly, scales as O(L). The thermodynamics and the ground 
state energy of the system have been obtained in this limit. Our exact solutions demonstrate 
explicitly the metal-insulator phase transition at half-filling, as the interaction of the local 
impurity moments and conduction band is switched on at zero temperature. At any nonzero 
temperature, the free energy has no singularity, indicating that there is no phase transition 
in this system. 
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